Total resummaion of leading logarithms vs standard description of the Polarized DIS 
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Total resummation of leading logarithms of x contributing to the spin-dependent structure func- 
tion gi ensures its steep rise at small x. DGLAP lacks such a resummation. Instead, the DGLAP ex- 
pressions for gi are complemented with special phenomenological fits for the initial parton densities. 
The singular factors x~ a in the fits mimic the resummation and also ensure the steep (power-like) 
rise of gi at the small- a; region. Furthermore, DGLAP by definition cannot describe the region of 
small Q 2 whereas our approach can do it. 

PACS numbers: 12.38.Cy 



I. INTRODUCTION 

The Standard Approach (SA) for theoretical investigation of DIS structure function gi(x, Q 2 ) involves DGLAPpJ 
and Standard fits|2| for the initial parton densities. In the SA framework, 

gDGLAP ig a 

convolution of the coefficient 

functions Cdglap and evolved (with respect to Q 2 ) parton distributions which are also expressed as a convolution 
of the splitting functions Pdglap and initial parton densities. The latter are found from experimental data at large 
x, x ~ 1 and Q 2 ~ 1 GeV 2 . As a result, SA accounts for the Q 2 -evolution through the DGLAP evolution equations 
whereas the x -evolution is accounted for through the fits which are found from phenomenological considerations. The 
reason for such asymmetric treating the Q 2 and x -evolutions in SA is that DGLAP was originally constructed for 
operating at large x where x- contributions from higher loops were small and could be neglected. In other words, the 
x -evolution can be neglected at large x. However, in the small- x region the situation looks opposite: logarithms of x 
are becoming quite sizable and should be accounted to all orders in a s . The total resummation of leading logarithms 
of x was first done in Refs. Q in the double-logarithmic approximation so that a s was kept fixed at an unknown scale 
and later in Refs. Q where the running a s effects were accounted for. Contrary to DGLAP where 

„DGLAP (i^fl\ (i\ 

a s = a s (Q ), (1) 

Ref. used the parametrization of a s suggested in Ref. [j| because the DGLAP parametrization of Eq. cannot 
be used at small x. The parametrization of Ref. 5] is universally good for both small x and large x. It converge to 
the DGLAP- parametrization at large x but differs from it at small x. 

Nevertheless, it is known that, despite DGLAP lacks the total resummation of In a;, it successfully operates at 
I < 1. As a result, the common opinion was formed that not only the total resummation of DL contributions in 
Refs. but also the much more accurate calculations performed in Refs. should be out of use at available x and 
might be of some importance in a distant future at extremely small x. In Ref. |(J we argued against such a point of 
view and explained why SA can be so successful at small x: in order to be able to describe the available experimental 
data, SA uses the singular fits of Refs. Q for the initialparton densities. Singular factors x~~ a in the fits mimic the 
total resummaton of Refs. Q. Using the results of Ref Q allows to simplify the rather sophisticated structure of the 
standard fits. 

Another essential difference between SA and our description of g\ is the obvious fact that DGLAP works at the 
kinematic regions of large Q 2 whereas our approach is valid for large and small Q 2 . The latter is important in 
particular for theoretical explanation of the COMPASS collaboration results. In Ref. we showed that g\ practically 
does not depend on x at small x, even at x C 1. Instead, it depends on the total invariant energy 2pq. Experimental 
investigation of this dependence is extremely interesting because according to our results g\, being positive at small 
2pq, can turn negative at greater values of this variable. The position of the turning point is sensitive to the ratio 
between the initial quark and gluon densities, so its experimental detection would enable to estimate this ratio. 
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II. DIFFERENCE BETWEEN DGLAP AND OUR APPROACH 



In DGLAP, gi is expressed through convolutions of the coefficient functions and evolved parton distributions. As 
convolutions look simpler in terms of integral transforms, it is convenient to represent g\ in the form of the Mellin 
integral. For example, the non-singlet component of g\ can be represented as follows: 



9lDGLAp{ x iQ 2 ) = ( e g/ 2 ) 



2m 



{I / xY C DGLAp{u)Sq{u) exp 



dk 



(2) 



with Cdglap{<jj) being the non-singlet coefficient functions, ~idglap{w, ct s ) the non-singlet anomalous dimensions 
and Sq(uj) the initial non-singlet quark densities in the Mellin (momentum) space. The expression for the singlet g\ 
is similar, though more involved. Both "/dglap and C'dglap are known in first two orders of the perturbative QCD. 
Technically, it is simpler to calculate them at integer values of to — n. In this case, the integrand of Eq. © is called 
the n-th momentum of g^ s . When the moments for different n are known, g NS at arbitrary values of oj is obtained 
with interpolation of the moments. Expressions for the initial quark densities are defined from phenomenological 
consideration, with fitting experimental data at x ~ 1. Eq. @ shows that jdglap govern the Q 2 - evolution whereas 
Cdglap evolve Sq(oj) in the a;-space from x ~ 1 into the smalkc region. When, at the x-space, the initial parton 
distributions Sq(x) are regular in x, i.e. do not — * oo when x — ► 0, the small- x asymptotics of g\ dglap is given by 
the well-known expression: 



N S S 
9l DGLAP > 9l DGLAP ~ ex P 



ln(l/x)\n(\n(Q 2 /p 2 )/Hp 2 /A 2 QCD ) 



(3) 



On the contrary, when the total resummation of the double- logarithms (DL) and single- logarithms of x is done 5], 
the Mellin representation for g^ s is 



du 



g? s (x, Q 2 ) = (e*/2) / —(l/xrC NS (fij)Sq(u) exp (H ns (oj) HQ 2 / p 2 )) , 



with new coefficient functions Cns, 



and anomalous dimensions Hns, 



Cns(oj) 



where 



H NS = (1/2) oj - V / ^ r 6R 



B(oj) = (4ttC f (1+oj/2)A(oj) + D{oj))/(2tt 2 ) 



(4) 



(5) 



(6) 



(7) 



D(u>) and A(oj) in Eq. are expressed in terms of p — ln(l/x), i] = \n(p 2 /K 2 QCD ), b = (33- 2n/)/127r and the color 
factors C F = 4/3, N = 3: 



2C F f°° _ up , ,p + TK 



(p + v) 



P + n 1 

— T — 



(8) 



A{oj) 



dpe 



Lrp + TT* J (p + -q) 2 +iT 
Hs and Cns account for DL and SL contributions to all orders in ot s . When x — » 0, 



9? S 



(x 2 /Q 



2\ Ajvs/2 



9? 



(x 2 /Q 



2\A S /2 



(9) 



(10) 



where the non-singlet and singlet intercepts are Ajvs = 0.42, As = 0.86. The x- behavior of Eq. (|10f> is much steeper 
than the one of Eq. © . Obviously, the total resummation of logarithms of x leads to the faster growth of gi when x 
decreasing compared to the one predicted by DGLAP, providing the input Sq in Eq. @ is a regular function of oj at 
oj -> 0. 
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III. STRUCTURE OF THE STANDARD DGLAP FITS 



Although there are different fits for 5q(x) in literature, all available fits include both regular and singular factors 
when x — > 0. For example, the typical expression is 

Sq(x) = Nr,x- a \{l - xf{\ + -yx 5 )] , (11) 

with N, rj being a normalization, a — 0.576, j3 — 2.67, 7 = 34.36 and <5 = 0.75. In the w -space Eq. (|llll is a sum of 
pole contributions: 

Sq(ui) = Nr] (w - a)- 1 + ^ m k {uj + Afc)" 1 , (12) 

with Afc > 0, so that the first term in Eq. Ill2[l corresponds to the singular factor x~ a of Eq. jH} . When the fit Eq. 1(11(1 
is substituted in Eq. the singular factor x~ a affects the small -x behavior of g\ and changes its asymptotics Eq. © 
for gi for the Regge asymptotics. Indeed, the small- x asymptotics is governed by the leading singularity w = a, so 



C(a)(l/ 2 ;) Q ((ln(Q 2 /A 2 ))/(ln( A1 2 /A 2 ))) 7(Q) . (13) 



9l DGLAP 



Obviously, the actual DGLAP asympotics of Eq. 1(13(1 is of the Regge type, it differs a lot from the conventional 
DGLAP asympotics of Eq. ® and looks similar to our asymptotics given by Eq. incorporating the singular 

factors into DGLAP fits ensures the steep rise of g DGLAP at small x and thereby leads to the success of DGLAP at 
small x. Ref. |(| demonstrates that without the singular factor x~ a in the fit of Eq. (|llfl . DGLAP would not be able 
to operate successfully at x < 0.05. In other words, the singular factors in DGLAP fits mimic the total resummation 
of logarithms of x of Eqs. 1(4110(1 . Although both 1)13(1 and 1)10(1 predict the Regge asymptotics for g\, there is a 
certain difference between them: Eq. 1(13(1 predicts that the intercept of gf s should be a = 0.57. As a is greater 
than the non-singlet intercept Ajvs = 0.42, the non-singlet g^ GLAP grows, when x — > 0, faster than our predictions. 
However, such a rise is too steep. It contradicts the results obtained in Refs. 4] and confirmed by several groups 
fitting HERMES data. Usually, the DGLAP equations for the non-singlets are written in the £-space as convolutions 
of splitting functions P qq with evolved parton distributions Aq and the latter are written as another convolution: 



Ag(a:) =C q (x,y)®8q{y), (14) 

with C q being the coefficient function. Written in this way, Aq is sometimes believed to be less singular than Sq 
because of the evolution. However applying the Mellin transform to Eq. 1(14(1 immediately disproves it. 



IV. gi AT SMALL x AND SMALL Q 2 

The COMPASS experiment measures the sing let gx at x ~ 10~ 3 and Q 2 < 1 GeV 2 , i.e. in the kinematic region 
where it is impossible to use DGLAP. Although formulae for singlet and non-singlet g\ are different, with formulae for 
the singlet being much more complicated, we can explain the essence of our approach, using Eq. for the non-singlet. 

In the COMPASS experiment Q 2 <C /i 2 . The expression for g\ at such small Q 2 with logarithmic accuracy is given 
by Eq. (0} where the Q 2 -dependence is dropped and x is replaced by fi 2 /2pq, so 



^—y/2pqrC« NS (cj)5q(u;) . (15) 

The expression for the gi singlet looks similar, though e 2 should be replaced by the averaged charge < e 2 > and 
CNs(u)$q((jj) should be replaced by the sum 

C g s (oj)Sq(u)+C 9 s (u)6g(u;) (16) 

so that 8q{io) and Sg(uj) are the initial quark and gluon densities respectively and Cg 9 are the singlet coefficient 
functions. Explicit expressions for Cg' 9 are given in Ref. 0. The standard fits for Sq and 5g contain singular factors 
~ x~ a which mimic the total resummation of leading logarithms of x. Such a resummation leads to the expressions 




FIG. 1: Gi evolution with decreasing z = /j 2 /2(pq) for different values of ratio r = 5g/8q: curve 1 - for r = 0, curve 2 - for 
r = —5 , curve 3 -for r = —8 and curve 4 -for r = —15. 



for the coefficient functions different from the DGLAP ones. After that the singular factors in the fits can be dropped 
and the initial parton densities can be approximated by constants: 

Sq ^N q Sg^N g , (17) 

so, one can write 

9i(Q 2 < M 2 ) /2)A g G 1 (z) (18) 

with 

/ lco dt > 
{l/zY[Cl + {N g /N q C g )] (19) 

where z = [i 2 /2pq. Obviously, G\ depends on the ratio N g /N q . The results for different values of the ratio r = N g /N q , 
Gi are plotted in Fig. ^ When the gluon density is neglected, i.e. N g — (curve 1), G\ being positive at x ~ 1, is 
getting negative very soon, at z < 0.5 and falls fast with decreasing z. When N g /N q — —5 (curve 2), Gi remains 
positive and not large until z ~ 10 _1 , turns negative at z ~ 0.03 and falls afterwards rapidly with decreasing z . This 
turning point where Gi changes its sign is very sensitive to the magnitude of the ratio r . For instance, at N g /N q = —8 
(curve 3), Gi passes through zero at z ~ 10~ 3 . When N g /N q < —10, Gi is positive at any experimentally reachable 
z (curve 4) . Therefore, the experimental measurement of the turning point would allow to draw conclusions on the 
interplay between the initial quark and gluon densities. 



V. CONCLUSION 



Comparison of Eqs. 10) and l|13l) shows explicitly that the singular factor x~ a in the Eq. 1|11[) for the initial quark 
density converts the exponential DGLAP-asympotics into the Regge one. On the other hand, comparison of Eqs. (|10|> 
and (|13|l demonstrates that the singular factors in the DGLAP fits mimic the total resummation of logarithms of 
x. These factors can be dropped when the total resummation of logarithms of x performed in Ref. Q| is taken into 
account. The remaining, regular x-terms of the DGLAP fits (the terms in squared brackets in Eq. (|llfl ') can obviously 
be simplified or even dropped at small x so that the rather complicated DGLAP fits can be replaced by constants. 
It immediately leads to an interesting conclusion: the DGLAP fits for Sq have been commonly believed to represent 
non-perturbative QCD effects but they actually mimic the contributions of the perturbative QCD, so the whole impact 
of the non-perturbative QCD on g\ at small x is not large and can be approximated by normalization constants. We 
have used the latter for studying the gi singlet at small Q 2 because this kinematic is presently investigated by the 
COMPASS collaboration. It turns out that g\ in this region depends on z — fi 2 /2pq only and practically does not 
depend on x. Numerical calculations show that the sign of g\ is positive at z close to 1 and can remain positive 
or become negative at smaller z, depending on the ratio between Sg and Sq. It is plotted in Fig. 1 for different 
values of Sg/Sq. Fig. 1 demonstrates that the position of the sign change point is sensitive to the ratio Sg/5q, so the 
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experimental measurement of this point would enable to estimate the impact of Sg. 
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